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Abstract 

The diffraction spectrum of the dart-rhombus random tiling of the plane is derived in rigorous terms. Using the 
theory of dimer models, it is shown that it consists of Bragg peaks and an absolutely continuous diffuse background, 
but no singular continuous component. The Bragg part is given explicitly. 
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1. Introduction 

Structure models of quasicrystals are usually based 
on the asumption of either energetical or entropi- 
cal stabilization of the material. Random tilings as 
possible structure models of the second kind where 
proposed || soon after the discovery of quasicrystals 
and studied thoroughly since, see e.g. jllj [l6| and 
references therein. Nonetheless, until today it is not 
yet clear which mechanism is dominating, although 
experiments are indicating a stochastic component in 
many cases fl3[| . Scaling arguments jlO|, [ll[] predict 
a singular continuous contribution to the diffraction 
spectrum for two-dimensional random tilings in ad- 
dition to the usual Bragg part and continuous back- 
ground. This should be visible in diffraction images 
of materials with so-called T-phases, see |lj and ref- 
erences therein, though it is not obvious how to dis- 
tinguish the different contributions. This underlines 
the necessity of investigating the diffraction of ran- 
dom tilings in more detail. 

In this article, we illustrate recently established re- 
sults |!| by the so-called dart-rhombus tiling. This 
two-dimensional model has crystallographic symme- 
tries and can be mapped onto the dimer model on the 
Fisher lattice. After introducing the tiling and the 
necessary mathematical tools, we calculate the two- 
point correlation functions and thereof the diffrac- 
tion spectrum. As in other crystallographic exam- 
ples, the spectrum can be shown to consist only of 
a Bragg part and an absolutely continuous back- 



ground, i.e. there is no singular continuous compo- 
nent. 

2. The dart-rhombus random tiling 

The dart-rhombus tiling is a filling of the plane, with- 
out gaps and overlaps, with 60°-rhombi of side I 
and darts made of two rhombus halves (Fig. [l]). In 




Fig. 1. The dart-rhombus tiling as dimer model on the 
Fisher lattice. The dots represent the atomic scatterers. 

addition to the usual face-to-face condition, we im- 
pose an alternation condition on the rhombi, such 
that neighbouring rhombi of equal orientation are 
excluded. Finally, to avoid pathological lines of alter- 
nating darts, we demand that two neighbouring darts 
must not share a short edge. These rules force the 
darts to form closed loops in a background of alter- 
nating rhombi. The minimal total rhombus density 
obviously is 1/3. This tiling can be mapped onto the 
fully packed dimer model on an Archimedian tiling 



known as Fisher's lattice in the context of statisti- 
cal mechanics. In order to control the densities of 
the different prototiles, we weigh them using activi- 
ties yi,Zi (Fig. ||), with Zi = e^ Mi etc., where [i i are 
chemical potentials and (3 the inverse temperature. 
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Fig. 2. Elementary cell of the Fisher lattice with activ- 
ities assigned to each bond and a typical random tiling 
with pi = 0.21, p 2 = 0.19, p a = 0.17. 

The grand-canonical partition function is there- 
fore equal to the dimer generating function. For any 
periodic graph with even number of sites in the ele- 
mentary cell, the latter can be computed as Pfaffian[| 
of the suitably activity- weighted adjacency matrix A 
. The calculation of the Pfaffian is simplified con- 
siderably by imposing periodic boundary conditions, 
but in the infinite volume limit the result holds also 
for free ones ^ Lemma 1]. 

Let us denote the rhombus densities by p. t (i = 
1, 2, 3) and the the dart densities by (j = 1, . . . , 6). 
There are several constraints on the densities. Closed 
dart loops require equal densities of opposite darts: 
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(1) 



Moreover, as each dart is accompanied by a corre- 
sponding rhombus, the remaining rhombi occur with 
equal frequency owing to the alternation condition, 



Pi ~ °\ = P2 ~ a 2 = P3 ~ CT 3- 



(2) 



Including the normalization constraint (the sum of 
the densities is 1), the number of independent param- 
eters (activities or densities) reduces to three. We 
exploit this freedom by setting all activities except 
Z\,Z2 and z% equal to 1. The dart-rhombus tiling 
undergoes second order phase transitions at 



1 + z\ + z\ + z\ = 2max{l, z\, zf , z%} or 

2,2,2 r> r 2 2 2i 

z i + z 2 + z 3 = ^maxjz-L, z 2 , z 3 } 



(3) 



with logarithmic (Onsager type) or square root di- 
vergence (Kasteleyn type), respectively. The point 
of maximum entropy is fixed by symmetry to p i = | , 
tjj = where darts and rhombi occupy half of the 
tiling area each. For further details see Jll| E^j. 

1 This is basically the square root of the determinant of an 
even antisymmetric matrix |15l Ch. IV. 2]. 



3. Diffraction theory 

For simplicity, we assume kinematic diffraction in 
the Fraunhofer picture B, i.e. diffraction at infinity 
from single-scattering. The diffracted intensity 7^ (a 
positive measure) is calculated as Fourier transform 
of the autocorrelation 7^ (see Q for details). It is 
known that every positive measure admits a unique 
decomposition into three parts /1 = p pp + /i sc + \i ac 
with respect to Lebesgue's measure, where pp, sc and 
ac stand for pure point, singular continuous and ab- 
solutely continuous p7| . In a diffraction spectrum, 
H pp are the Bragg peaks and /i ac the usual diffuse 
background or Laue scattering. A singular contin- 
uous part can be encountered in ID substitutional 
sequences, cf ||, and is also expected for 2D qua- 
sicrystalline random tilings |l0|, [n], E|. 

Consider the so-called weighted Dirac comb [|| 



w{x)5 x 



(4) 



on a lattice f , where S x is the unit point mea- 
sure (Dirac measure) concentrated at x, and w{x) G 
{0, 1} is chosen in order to obtain a specific member 
of the random tiling ensemble. Its autocorrelation 
7c; is (almost surely) 



7a 



v(z)5 z 



(5) 



Here, A = T — r = Tis the set of difference vec- 
tors and the autocorrelation coefficient v{z) can be 
calculated according to 



v(z) = lim 



1 



-fi^oo voI(-Br) 



^2 w(y)w(y + z), (6) 



where Br is the ball of radius R around the origin, 
Ajj = Ad Br and — the complex conjugation. Thus, 
v(z) is simply the probability of having two scatter- 
ers at distance z, which a.s. exists. 

We use standard Fourier theory of tempered dis- 
tributions, for our conventions see 

4. Diffraction of the dart-rhombus tiling 

We decorate the tiling with point scatterers 5 X ac- 
cording to Fig. 0. More realistic atomic profiles can 
be handled with the convolution theorem [[l7|, Ch. 
IX] . The scatterers may have complex strengths h Pi , 
resp. h G . , where the strengths of opposite darts are 
supposed to be equal. This constraint simplifies the 
calculations but is not necessary. The point set of 
all possible atomic positions is a Kagome grid with 
minimal vertex distance 1/2. We write it as tri- 
angular lattice r with a rhombic elementary cell E 
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containing the nine scatterers positions for the dif- 
ferent tiles (cf Fig. ||). Introducing basis vectors 
ei = (-v/3,0)*, e 2 = 1/2(^3,3)* for T and E, the po- 
sitions p in E (with corresponding density) are given 
by (o = 1/4(V5, 1)*, b = 1/4(^3,-1)*) ' 



p Pi =3a, 

=a, 
Pa d = 3a 



•6, 



Pp 2 = 2a - 6, 
p CT2 = 2a + 6, 
Per, =5a, 



Pp 3 = a + b, 

p„ 3 =3a-b, (7) 

p a = 4a-b. 



We now have to calculate the autocorrelation or 
the joint occupation probability of the dimers. As 
we will see, the autocorrelation coefficients can be 
split into a constant term and one decreasing with 
the distance between the scatterers. After taking 
the Fourier transform, the first will yield the Bragg 
peaks whereas the second will be responsible for the 
continuous part of the spectrum; we will show that 
this can be represented by a continuous function and 
hence contains no singular contribution. 

Using Gibbs' weak phase rule one can prove 

the ergodicity of the model and thus justify the cal- 
culation of the diffraction spectrum via the ensemble 
average. 

Let rj kk , be the occupation variable that takes the 
value 1 if the bond between k and k' is occupied and 
otherwise; fj kk , := 1 — rj kk ,. As was shown in 
the probability P a p of bonds a and being occupied 
simultaneously is given by 



P, 



a/3 



+ (Vk, 



)( r lk l3 k' l3 ) 



(Vk n 



(8) 



P a Pfi 
-A 



i 



), 



with p a the density of dimers that can occupy the 
bond (k a k' a ) and A the weighted adjacency matrix. 

We consider the constant part of the autocorrela- 
tion first. Combining (|^) and (||) we get for the point 
set of scatterers with density 2 / V3 
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where * denotes convolution. Computing the Fourier 
transform with Poisson's summation formula ^ Eq. 
12] using (Q) and (|l|), the pure point part of the spec- 
trum is (almost surely) 



(7), 



E 

(M)er* 



KiPl + (~ 1 ) kh P2P2 + (-l)\ 3 />3 



2 cos 



, Tv(k+l) 
3 



\k+li 



(9) 



{-l) l K 2 a 2 + {-l) k K 3 u^j 



where T* is spanned by e\ = (-^, — |), e 2 = (0, |). 

It remains to calculate the other part of @ . Ay. y 
is nonvanishing only if k a and k' a are connected; in 
this case A^. k' = ez a , with e = ±1 according to the 
direction of the arrows in Fig. ||. 

Since A is the adjacency matrix of a graph that 
is an (m, n)-periodic array of elementary cells with 
toroidal boundary conditions and therefore cyclic, it 
can be reduced to the diagonal form A = diag{Aj} by 
a Fourier-type similarity transformation with matrix 
Skk' = {mn)^ 1 / 2 cxp(2Tri(k 1 k' 1 /m + k 2 k' 2 /n)). A^ 1 
is then determined by B 



(m,n) 

E s^xfslj. (10) 

3=(1,1) 



In the infinite volume limit, the sums approach inte- 
grals (Weyl's Lemma), and by introducing r = k'—k, 
(p 1 = 2 r Kij 1 /m etc. we obtain 



A" 1 = J- 

^kk' 47r 2 



A-^^y-oV^. (11) 



To determine A , observe that the inverse of 



A = 



can be computed easily by any computer-algebra 
package but unfortunately does not fit onto this page. 

Defining the coupling function [x, y] PlP2 for two 
dimers in elementary cells at distance r = xei -\-yei, 
with {x,y) £ Z 2 , where the dimers occupy positions 
pi, resp. p2 in each elementary cell, we rewrite ( fli"| ) 
in more explicit form 
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det (X(ip 1 ,ip 2 )) 

where the determinant of A is given by 

det (A) = a + 2b cos ip 1 + 2c cos ip 2 + 2d cos^-l — ip 2 ) 



with a = z\ + z| + z| + l, b = z\z\ — zf, c = z\z\ — z\ 
and d = z 2 z 2 — zf; g can be taken from the corre- 
sponding entry in A -1 and is finite. In order to de- 
termine the spectral type of the diffraction, we are 
interested in the asymptotic behaviour of [x, y] PlP2 
for large r . Substituting v = e~* v i and w — e~ lV2 , 
we obtain that / = 47r 2 [x, y) PlP2 is 

—g(pi,p 2 ,v,w)v~ x w~ v dvdw 

v 2 (bw + d) + v(aw + c(w 2 + 1)) + w(b + dw) ' 

S 1 xS 1 
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displays Bragg peaks only. For the Onsager case, we 
re-substitute w = e~ lip2 in (|l3|). Because of the sym- 
metry of the kernel, it is sufficient to integrate from 
to 7r. The kernel reaches its maximum value 1 only 
at ip2 — or 7r and remains smaller elsewhere. Us- 
ing the same argument as in the treatment of the 
lozenge tiling in [||, we estimate the kernel by a de- 
creasing/increasing straight line. From the result- 
ing asymptotic behaviour we conclude that the dif- 
fuse part of the spectrum is an absolutely continuous 
measure as well. 



Fig. 3. Diffraction image of the tiling of Fig. Q with scat- 
terers of equal strength 1. It was calculated numerically 
by means of standard FFT, because this is simpler than 
using the exact expression for the ac part. 



with S 1 the unit circle. We integrate over v for x < 0; 
the case of positive x can be treated analogously. The 
integrand is singular at 



v± 



-a±sja 2 -Aj3 
2(bw + d) 



(12) 



with a — (aw+c(w 2 + 1)) and (3 = w(bw + d)(b+dw), 
but only v+ lies inside the unit circle. Thus, 



/ = 2m 



S 1 



g(p 1 ,p 2 ,v+,w)v + x w y 
(v + — V-)(bw + d) 



dw. 



(13) 



Away from the phase transitions, it can be shown 
that \v + \ < 1. With v + = max ¥ , 2 w+ we get 



III < 2tt 



\g(p 1 ,p 2 ,v+,w)\\v + \ 



-dw 
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< 27r|6+|" 
= O (e-* 1 



v + — v-)(bw + d)\ 

\g(p 1 ,p 2 ,v+,w)\ 
si \(v+ - v-)(bw + d)\' 

*n (u) 



-dw 



for some positive constant t\ , because the remaining 
integral stays finite. Since the coupling function is 
invariant under interchange of x and y, this can be 
shown for y as well. As I is maximal for x = for 
arbitrary but fixed y and vice versa, we conclude that 



[ x : y]pu 



O e 



-(ti\x\+t 2 \y\) 



(15) 



The non-constant part of the correlation function 
in (||) consists basically of products of [x,y]. With 
such an asymptotic behaviour, one can show that its 
Fourier transform indeed converges towards a con- 
tinuous function on R 2 /r (cf ||, Addendum]). 

At the Kasteleyn phase transitions, we get crys- 
tals consisting only of one rhombus orientation and 
the corresponding darts (cf |10] ) . The spectrum thus 
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